Vortex dynamics of a d + is- wave superconductor 
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The vortex dynamics of a d + is-wave superconductor is studied numerically by simulating the 
time-dependent Ginzburg-Landau equations. The critical fields, the free flux flow, and the flux flow 
in the presence of twin-boundaries are discussed. The relaxation rate of the order parameter turns 
out to play an important role in the flux flow. We also address briefly the intrinsic Hall effect in d- 
and d + is-wave superconductors. 
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In recent years the symmetry of the pairing function 
has been one of the interesting topics in the field of high 
temperature superconductor. It is widely accepted that 
the dominant pairing wave function has a d x i_ y i symme- 
try, as supported by experiments using phase-sensitive 
devices, such as Josephson junctions or SupercjQnduct- 
ing Quantum Interference Devices (SQUID). uu How- 
ever, the subdominant pairing channels |Sik1i as s-wave 
or d xy -w&ve channels are still possible. Era The mixed 
state of s + d superconductors was first discussed by 
Ruckenstein et—al. u and_that of d + is superconduc- 
tors by KotliarEI In Ref.|ij], it is pointed out that the 
resonating-valence-bond mechanism can lead to s-wave 
and d- wave- like Cooper pairings, and a mixture of s and d 
waves with a well-defined relative phase close to 9 = it/2 
is energetically favored. Importantly, the superconduct- 
ing state is time-reversal-symmetry (T hereafter) break- 
ing and the energy gap is nodeless unless 9 = or n. 
Moreover, surfaces and interfaces, grain boundaries, and 
other pair-breaking defects haip,all been shown to en- 
hance the T-breaking states. EJE 2 ] 

The vortex structure and vortex dynamics of d-wave 
superconductor haV|e-,becn studied numerically in detail 
by several methodai3ilZI. Previous simulations for the 
d+is- wave superconductors showed that the vortex struc- 
ture of d + is- wave superconductor is different from that 
of the d-wave superconductorEa: The spatial profile of 
the moduli of s- and d-wave components in the d + is- 
wave state has a two-fold symmetry, in contrast to the 
four- fold symmetry of the magnetic-field-induced s-wave 
component of d-wave superconductor. However, it turns 
out that such a d + is state can only be stabilized at 
extremely low temperatures. With increasing tempera- 
ture, the amplitude of s-wave component decreases and 
its symmetry changes from two-fold to four-fold. Here we 
would like to extend our study to the vortex dynamics of 
d + is superconductors. For conventional superconduc- 
tors the free-flux-flow (FFF) resistance is linear in the 
magnetic field induction B (with B <C B C 2). The situa- 
tion is not clear yet for high-T c and other unconventional 
superconductors. Because of the multiple components 
of the order parameter, this topic is highly nontrivial. 



We also study the vortex motion in the presence of a 
twin-boundary and the intrinsic Hall effect of d + zs-wave 
superconductors. 

We start with a model of an isotropic two-dimensional 
Fermi liquid with attractive interactions in both s- and 
d- channels. Obviously, when only one of the two inter- 
actions is attractive, the ground state is a pure state with 
the appropriate pairing symmetry. When both channels 
are attractive, the competition will lead to either a pure 
or a mixed pairing function. The GL theory for a su- 
perconductor with two attractive channels has been pre- 
sented by Ren, Xu, and TingtJ based on Gor'kov equa- 
tions. Assuming pure dissipative dynamics, the GL equa- 
tions can be expressed as follows: 
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In these equations, the two order parameters, S and 
D, are normalized by A = y/A/Zahii/Td/T) with a w 
7£(3)/8(7rT c ) 2 , the space by the coherence length £, and 
the vector potential A by &o/2irt; with $o = h/2e be- 
ing the flux quantum, respectively. In Eq.IlL-p^ may be 
expressed as a function of temperature Tiij'EII 



a s = lii(T s /T)/\n{T d /T) 



(4) 



where T s and T d may be viewed as the apparent su- 
perconducting transition temperatures for the s-wave 

i 

and the d-wave respectively, with T s oc e N <°)^ and 



Td (x e N< -°> v d . ad = 1 unless specified otherwise. Here 
N(0) is the density of states at the Fermi surface, V s 
and Vd are the effective attractive interaction strengths 
in the s- and d-wave channels, respectively. II = iV + A, 
life = ifcllfe. k is the GL parameter. The time t is nor- 
malized by t = cr„£ 2 with a n the normal-state conductiv- 
ity of the superconductor. rj St d — T s _d/r, with T s ,d being 
the relaxation time of the s- and d- wave order parame- 
ters, respectively. For gapless conventional superconduc- 
tors, the relaxation time of .the order parameter near T c 
is t s = Trh/8k(T c — T) oc ^ 2 EJ For unconventional super- 
conductors, it is plausible to assume the general scaling 
law T s> d ~ £, z with z as the dynamical exponent, which 
should hold at least near the critical temperature. In 
this case, it is clear that r] St d ~ £ 2 " 2 would depend on 
temperature if z ^ 2. However, in the absence of a mi- 
croscopic theory on the relaxation times, we shall content 
ourselves to treat them as phenomenological parameters, 
the effect of which on the vortex dynamics is examined 
in this paper. 

The scale for the magnetic field is set by the upper 
critical field H C 2 . Because of its nontrivial nature in the 
present system, we would like to discuss it before moving 
to the dynamical properties of the vortices. Following the 
same method described by Sigrist et a/E3, and using the 
relation [n x ,n a ] = iB 1 we introduce a pair of operators 
a = -^=(U X + illy), a + = -jL=(Il x - ill y ), and have the 

commutative relation [a, a + ] = 1. Substituting a,a + into 
the linearized Eq.(|l|) and Eq.(||) in the static case, we 
obtain: 



(-2a s + 2B(l + 2h))S + B(aa + a+a + )D = 0, 
(-1 + B{1 + 2n)D + B(aa + a + a+)S = 0, 



(5) 



£ 



suming W 



where n = a + a. Therefore, the order parameter can 
be expanded in this occupation representation as $ = 

, m v I, where n I m) = m \ m). Eq.da) can 
m=0 \ d ™ I m)d J [ ^ 

not be solved exactly, so we treat it variationally by as- 
so | 0) s \ / S2 I 2). \ ( S4 I 4) s ■ 
do | 0> d J + \d 2 | 2) d ) + \di | A) d 
H C 2 is determined by the condition that the ground 
state energy of the eigenvalue problem Eq.(|^) is zero. 
It turns out that the variational ground state is W = 
W S2 |2) S W 
d \0) d ) + { J + ^ d A I 4> d 
sponding H C 2 is given by the root of the equation 



, and the corre- 



60B 3 - 86B 2 + 9a s B 2 + (10 - 10a s )B + a s = 0. (6) 

The temperature dependence of H C 2 calculated in this 
way is shown as the solid line in Fig.(l). Here we use 
Eq. (0) for the temperature dependence of a s , and assume 
T d = 100K and T s = 90K. In order to check the relia- 
bility of the variational treatment, we also simulate the 
upper critical field of such a d+is superconductor numer- 
ically by solving Eq.ffl) ~ Eq.(0) with the finite element 
method. In practice, we fix the magnetic field induction 



B by specifying one or two vortices in a square unit cell 
with periodic boundary conditions. The side length is 
varied so as to change the magnetic induction B. The 
technics] details of the simulation have been given else 
whereO The external magnetic field H can be derived 
from the Virial relation.E3 With H the Gibbs free energy 
density can be constructed, as shown in Fig. 2. We can 
read off H c \ from the intersection of Meissner state line 
and the mixed-state line, and H C 2 from the intersection 
of the mixed-state line and the normal-state line. The 
simulation result for the upper critical field is shown as 
black squares in Fig.l. The variational result (solid line) 
appears to be in excellent agreement with the simulation 
result, indicating that expanding the trial wave function 
up tOr-the occupation state of |4) is already rather re- 
liable £-3 Note that H C 2 is larger than that of the pure 
d-wave superconductor, which is always Bq = $o/27r£ 2 
in our case. From Fig. 2, it is obvious that the magnetic 
filed-induced transitions at H — H c \ and H = H C 2 are 
both of the usual second order transition, which should 
be compared to the unusual first order transitions found 
numerically in some p-wave superconductors EB 
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FIG. 1. The upper critical field H C 2/Bq as a function of 
temperature T with Bo as the upper critical field of pure 
d-wave. The solid line is the approximation result determined 
from Eq.(6). The scatter point is the numerical result. Here 
T s is set to 907^, T d is set to 100 K. 
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FIG. 2. Free vortex flow resistivity as a function of the 
applied field for r? s = 1, 10, and 0.01. Here a 3 = 0.97, 0.85, 
0.67, and —1 in (a), (b), (c), and (d) respectively. 

We now discuss the vortex flow driven by a transport 
current J in the mixed state (i.e., with B < H C 2). This is 
realized by by requiring V x H e = J in Eq.(J3|). Here we 
have chosen a gauge in which the electrostatic potential 
does not appear, so that the local electric field is simply 
E = —d t A. To investigate the relaxation effects of the 
order parameter on the FFF resistivity, we chose three 
values of the relaxation coefficient of the order parameter 
S and D as :(1) rj s = 2 m = 1, (2)r) B = 2% = 10 (> 1), 
(3) r] s = 2rj d = 0.01 (<C 1). The field dependence of the 
FFF resistivity is shown in Fig. 3. Noticeably it evolves 
from a convex to a concave with increasing ?y s , albeit 
slight difference exists at the four temperatures shown 
in the four panels. This is the correct trend by general 
reasoning: The motion of vortices is [equivalent to the 
phase slipping of the order parameter ,E3 so that a small 
relaxation time means a quick rate of phase slipping, and 
thus a large resistance. Thus, with increasing magnetic 
field the vortex flow resistivity approaches the normal 
state resistivity more quickly. It follows from Fig. 3 that 
the effect of the relaxation time is more prominent at 
lower temperatures (or larger a s ). 
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FIG. 3. Free vortex flow resistivity as a function of the 
applied field for r? s = 1, 10, and 0.01. Here a s = 0.97, 0.85, 
0.67, and —1 in (a), (b), (c), and (d) respectively. 

Next, we investigate the vortex motion of the d + is su- 
perconductor in the presence of twin-boundaries and then 
look into the pinning effect. A periodic array of twin- 
boundaries (with a transverse spacing of L — 10. 8£) are 
assumed and described by a, = ai,o — u i Yl 3{y~yo — kL), 

k 

where the subscript i stands for s or d. Here, u% describes 
the variation of a St d across the twin boundary along the 
line y = yo due to the local mis-orientation or chemi- 
cal contamination. Along the twin-boundary, we apply 
a transport current J = Jx. The vortex motion will be 
pinned by the twin-boundary up to a depinning current 
J = J c , which we shall determine. In the following sim- 
ulation we fix rj s — 2r\d = 1. Figure 4 gives the current 
dependence of the flux- flow resistivity, which turns out to 
be highly nonlinear. In the present simulation, we may 
expect a simple result of the overdamped vortex motion 
in a periodic pinning potential: pj Pnr=j a-\/l — (J c /J) 2 
at J > J c (the solid lines in Fig. 4), 113 where a is the 
asymptotic reduced resistivity and J c can be thus deter- 
mined. We see that a higher depinning current arises 
from higher values of u s (or Ud)- This is because of the 
increasing suppression of the the amplitudes of s- wave 
and d- wave components at the boundary. In fact, with 
a suppression of the order parameter at the twin, the 
vortices loose less condensation energy, and thus have a 
lower energy than they would do in the bulk. For a d- 
wave superconductor, u s is irrelevant because a s < 0, 
only Ud takes effect, so it is reasonable that a d- wave su- 
perconductor has a lower depinning current than a d + is 
superconductor does with the same twin boundary. 
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FIG. 4. The current dependence of the flux-flow resistivity 
in the presence of a twin-boundary. The solid lines represent 
the fitting at low J(see the text). 



Finally, we turn our attention to the so called intrinsic 
Hall effect, which occurs for anisotropic vortex structures 
even without considering normal Hall conductivity a 7 ^. 
Alvarez et al. have found that in d- wave superconductors 



the intrinsic Hall effect depends on the orientation angle, 
if, of the driving current taughly as ~ sin(&<p), and in- 
creases nonlinear ly with J £3 It is believed that this effect 
is due to the four-lobe structure of the d- wave vortices. 
So here we would like to explore the effect of the two- 
fold symmetry of s- and d- components on the intrinsic 
Hall effect. As usual, the intrinsic Hall effect is mea- 
sured by tan## = Eav,-l/Eav,\\i where E^y and Eav,± 
are the components of E^y perpendicular and parallel 
to the applied current, respectively. Figure 5 presents 
tan0# versus tp for (a) a d- wave superconductor with 
a s = — 1.0 and ay — 1.0; (b) a d + is- wave superconduc- 
tor with a s = 0.85 and a d = 1.0. Here B = 0.034B c2 , 
and |J| = 0.08 at which p/p n ~ 0.6 in both cases (a) and 
(b). It can be seen that for the d- wave superconduc- 
tor (Fig. 5(a)), tan## has two peaks (at ip p = 15° and 
ifp = 75°), and is identically zero at ip = 0°,45°, and 
90°. For the d + is superconductor, tan6*_y is still zero at 
ip = 0°,45° and 90°, but there are new sign reversals at 
approximately ip — 22.5° and 67.5°. Clearly, the vanish- 
ing of tan Ofj at (/? = 0°,45° and 90° can be attributed 
to the fact that the vortex(and thus the supercurrent 
around it) can adjust its symmetry axis to the the direc- 
tion of the driving force from the the applied current. At 
other directions of the driving force, the vortex can only 
partially adjust its symmetry axis. The new sign rever- 
sal for the d + is superconductors is more subtle. In this 
case, the vortex profile is two-fold symmetric. Therefore, 
there are two non-equivalent configurations for the vor- 
tex to have its (reflection) symmetry axis along the x, y, 
or the diagonal directions. We suspect that the abnormal 
sign reversal results from a switching between these two 
meta-stable configurations as the direction of the driving 
force changes. 



:/v 


■---« 




(a) 


^ 




-■ ■-.. 


w 


f f\ 


'•• 


-. f 


(b) 


•• * 


J 

• 


'••I 


\i 



FIG. 5. tan 6h as function of the current orientation angle 
ip at |J| = 0.08 and B = 0.034B c2 . (a) a s = -1.0 and 
ad — 1.0; (b) a s — 0.85, ad = 1.0. The lines are guides for 
the eye. 

In summary, the dynamics of the vortices in both d+is- 
wave and d-wave superconductors is studied. The upper 
critical field of the d+is- wave superconductors is studied 



analytically and numerically. From simulation results, 
the curvature of the FFF resistivity as a function of the 
magnetic field strongly depends on the relaxation rate 
of the order parameter. The flux flow in the presence 
of a twin boundary is also addressed. Finally, the in- 
trinsic Hall effect of d- and d + is- wave superconductors 
are studied numerically . We find the orientation depen- 
dence of the intrinsic Hall effect in these two types of 
superconductors are rather different. 
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